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Search for the Lepton Number Violating Decay r ^ 

Chris Brown** for the BABAR collaboration 

'^Department of Physics and Astronomy, University of Victoria, 
P.O. Box 3055, STN CSC, Victoria BC V8W 3P6, Canada 

Using data collected with the BABAR detector between 1999 and 2001, we describe a preliminary search for 
the neutrinoless decay t —> fi'y. This data sample includes data collected both on and off the T(4S) resonance and 
corresponds to 56 x 10^ produced tau pair events. The search has an efficiency of 5.2 ± 0.1{MCstat) ± 0.5{sys)% 
and an expected background rate of 7.8 ± 1.4 events. We select 13 events in the final sample. As there is no 
evidence for a signal in this data, we set a preliminary upper limit of B(r fj.'y) < 2.0 x 10~®@90% CL. 



1. Introduction 

The decay t — > ^7 is an anticipated 
lepton-number violating process in supersym- 
metric models |^,^,^, left-right supersymmetric 
models [1 and in supersymmetric string unified 
models[5 . For some ranges of model parameters, 
decay rates as high several parts per million are 
expected for this decay |^,||, even in light of the 
current experimental limit on the related /i — > 67 
decay Assuming that neutrinos have mass, the 
standard model branching ratio is expected at the 
0(10^'^'') level. The current best published limit 
is B(r ^7) < 1.1 X 10-s from the CLEO ex- 
periment using 12.6 million t pairs Q. 

This analysis uses 56 x 10^ r+r^ events from 
the 1999-2001 BABAR data set to give a prelim- 
inary result of B(r — > fij). The general signature 
of the T — > /X7 signal is the presence of an iso- 
lated ^ and 7 which have an invariant mass con- 
sistent with that of the t( 1.777 GeV) and with an 
energy in the centre-of-mass system (CMS) con- 
sistent with the beam energy in the CMS (5.29 
GeV) and the rest of the particles in the event 
having properties that are consistent with being 
produced in a generic 1-prong t decay^. The anal- 
ysis is performed using a blinded mass-energy re- 
gion in the data which corresponds to a 3cr error 
ellipse centred on the expected peak of the distri- 
bution. 



'1-prong' decay contains a single charged particle 
amongst its decay products. 



A cut-based approach is used, in which the 
background in the signal-box is estimated from 
extrapolations using side-band data and verified 
with the Monte Carlo simulation. The determi- 
nation of selection criteria in this analysis was 
based on Monte Carlo simulation of background 
and signal as well as on data in side-band regions 
further away from the signal-box. 

Background sources arising from non-r+r^ 
sources, the most problematic of which are radia- 
tive fJ,'^ ^jL~ pair events, are reduced to small levels. 
This leaves the generic 1-prong r decays as the 
major source of background. The dominant, and 
irreducible, background is from radiative muonic 
decays of the t in which the two neutrinos have 
very little energy. 

2. Monte Carlo Samples 

The Monte Carlo simulation uses a complete 
description of the BABAR detector response em- 
ploying the GEANT4§| software. The most im- 
portant simulation samples are those involving 
e~^e~ T^T~ decays. These include 27.6 million 
generic r+r^ events, which use KORALB em- 
ploying the TAUOLA decay package]^, and the 
r — > /17 signal incorporated into TAUOLA. Forty 
thousand signal events were generated which form 
the bulk of the signal Monte Carlo used for the 
analysis. A second signal sample with differing 
beam background contributions was produced in 
order to provide a cross check on the sensitivity of 
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the selection to different background conditions. 

Additionally, simulated z^"*"^" events, which 
uses the AfkQed generator ||l0|, provide tools 
for ^ particle identification and some guid- 
ance on e^e^ l-J'^ 1^^ event suppression. 
Other potential backgrounds include e+e" — > 
uu; dd; ss, e~^e~ — > cc and e^e~ bb which were 
studied using the corresponding generic BABAR 
simulated event samples. As the selection is re- 
stricted to events in which both t's decay via 1- 
prong modes, the e+e~ — > bb background is com- 
pletely negligible as indicated by the simulated 
events. The e+e~ — > cc and light-quark back- 
ground contributions are also found to be very 
small, with no events in the simulation surviv- 
ing from these two sources after all cuts have 
been applied, and consequently have not been in- 
cluded further in the data-Monte Carlo compar- 
isons. The two-photon processes are negligible 
backgrounds for events in which the beam energy 
is detected in a CMS hemisphere. 

3. Mass and Energy Determination 

Because there is no missing energy in the r — > 
/j,7 signal, the mass of the measured decay prod- 
ucts is the r mass and the energy is the full energy 
of the r. Ignoring the effects of initial-state ra- 
diation, the energy of the final state particles is 
equal to the full energy of the beam in the CMS. 
This situation lends itself very well to the use of a 
mass of the fi'y system calculated from the kine- 
matic fit employing the beam energy constraint 
on the energy of the system as used by ARGUS, 
and denoted as iriEc- The energy variable used 
is the difference between the measured fij energy 
and the beam energy, denoted by the symbol AE. 
The distributions of these variables are shown in 
Figure ||(a)-(c). The resolutions of the core of 
these distributions, which represent those events 
with well reconstructed photons and tracks, are 
(TAE = 88 MeV and am^c = 19 MeV. 

However, as is evident by the diagonal band of 
events present in Figure |^(c), when initial state 
radiation shifts the true energy of the r from the 
beam energy, there a negative correlation between 
AE and niEc and a substantial loss of resolu- 
tion. The low-energy tail of the AE distribu- 
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Figure 1. Distributions of (a) AE (b) niEC and 
(c) TriEc vs AE for the Monte Carlo Simulation 
of the r ii"f signal. 



tion for well reconstructed ttiec j evident in Fig- 
ure |l|(c) , is populated by events with a photon re- 
constructed with the correct direction, but with 
significant energy loss. 

Figure ^ shows the various regions in the rriEC 
vs AE plane used for determining the selection 
criteria, measuring the background, and defin- 
ing the signal and blinded regions. A 'Grand 
Side-Band' is defined as the region within the 
1.0 <AE< 0.5 GeV and 1.5 < ttiec < 2.1 GeV 
bounds. It is used for evaluating the reliability of 
the estimation of the signal efficiency and this is 
also shown on the figure. 

The signal-box is an elliptical region centred 
on the peak of the two-dimensional distribution 
as determined by the Monte Carlo. For the posi- 
tive side of AE, the ellipse has a 3cr half-axis for 
both the AE and rriEC axes whereas for the neg- 
ative side of AE the ellipse has a 2a half-axis in 
AE and a 3cr half- axis for niEC- The resolutions 
used in defining the signal-box are those of the 
core Gaussians obtained from the signal Monte 
Carlo. This asymmetric shape provides an op- 
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Figure 2. Regions of the niEC vs AE plane shown 
schematically overlaying the distribution of the 
simulated signal. The signal region is the assy- 
metric ellipse within the elliptical blinded region. 



timal signal-box given the resolution in AE and 
the presence of increasing background in the neg- 
ative AE region that is not present for positive 
AE. 

Events that are ultimately selected and in- 
cluded in a study of the niEC and AE variables 
must first pass a number of selection criteria. 
These requirements are designed primarily to re- 
duce the number of Bhabha and n'^ n~ events. 
The requirements of the standard BABAR back- 
ground filter which the signal events pass: 

• 2 charged tracks 

• zero net charge 

• bi| + \P2\ < 9.0 GeV 

• Ei + E2< 5.0 GeV 

• Ei/\pi \ < 0.8 OR E2/\p2\ < 0.8 

• EcM - (|pi|cA/ + \p2\cm) > 

• Pt{CM)/{EcM - {\pi\cM + \p2\cm)) > 0.07 

• charged track separation > 90° in the CMS 



• at least one 7 

• fiducial region: 0.775 < cos Otrack < 0.940 

where Pi(2) is the track momentum; -Ei(2) is 
the calorimeter energy associated with the track; 
Pt{CM) is the total transverse momentum of 
charged tracks in the CMS; and subscript CM 
indicates when the quantity is in the CMS. 

These requirements have an efficiency of 37.3 ± 
0.2%. These efficiency losses arise largely from 
the fiducial acceptance of the detector, but also 
from the branching ratio to 1-prong modes. At 
this stage of the analysis, 22 x 10^ events are se- 
lected in the data. Monte Carlo studies indicate 
that 83% of the signal events pass these require- 
ments for events which would have otherwise been 
selected. 

Subsequent selections employ requirements to 
remove the non-r background sources and the 
e^e^ —> fi'^fi~^ and e^e~ —>■ e^e~ events in 
particular. Cuts include a specific requirement 
on the observed CMS energy, Elfj2E^^'^^ < 
0.95, to exclude the aforementioned background 
sources in general and a /i veto for reduction of 
di-muon backgrounds in particular. The event is 
tagged via an electron or ft, > Ivr*' on the oppos- 
ing side. The missing mass on the tag side is 
restricted to reduce backgrounds from lost or dif- 
ficlt to reconstruct tracks. To be consistent with a 
decay through an intermediate r, the tag track is 
required to have less than 80% of the CMS beam 
energy. In general, this will not be the case for 
the non-r backgrounds. 

On the signal side, we impose the requirement 
of an identified /i and an associated 7 with energy 
in excess of 400 MeV. These cuts are applied in 
addition to the obvious m ec and AE cuts for the 
final selection. 

The cut progression shows that once the non-r 
background is removed by requiring an electron- 
VS-/U7 or ft > l7r°-vs-/i7 events, the r Monte Carlo 
simulation tracks the data reasonably well. 

We present in Figures ^ and ^ a sample of the 
distributions of variables employed in the analysis 
after all other requirements, apart from that us- 
ing the one shown, have been applied. The data 
distributions are represented by points and the 
r Monte Carlo simulation by histograms in these 
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figures. The normalizations of the Monte Carlo 
distributions is fixed by the luminosity. 
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Figure 3. Distributions of niopp after all other 

requirements, apart from the final rriEC-^E, have 
been applied. The data and t Monte Carlo are in 
the top plot and the t ^ ^'y signal Monte Carlo 
in the bottom plot. The parts of the distributions 
accepted in the selection are indicated. 



These distributions arc well described by the 
T~^T~ simulation in both shape and overall nor- 
malization. There are 604 events in the data in 
these distributions and the ratio of data to the 
Monte Carlo expectation is 1.022 ± 0M9{stat) ± 
0.025(norm). The first error is statistical and 
the second is the systematic error associated with 
normalization of the Monte Carlo events. This 
normalization uncertainty is a component of the 
systematic error associated with the efficiency of 
selecting the signal. This includes uncertainties 
in integrated luminosity, cross-section, radiation 
treatment in the generator and branching ratio 
uncertainties. 

Because the signal has a relatively stiff momen- 
tum spectrum, muon particle identification (/x- 
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Figure 4. Distributions of E-^ after all other re- 
quirements, apart from the final itlec-^E, have 
been applied. The data and t Monte Carlo are in 
the top plot and the t ^ /j.'y signal Monte Carlo 
in the bottom plot. 



PID) efficiency using a sample of radiative n pair 

events {e^e~ — > fX^ jjT^) from the same data set 
used for the rest of the analysis was exploited 
to study the /x-PID efficiency. Contamination 
from sources other than the /i pair, such as t~^t~ 
events, are negligible in this sample. 

The ratio of the efficiencies as a function of 
LAB momentum of the data to that of the 
e+e" — > simulation was studied and a 

correction based upon the control sample was ex- 
tracted. Similarly obtained is the correction for 
the T+T~ events in the Monte Carlo sample and 
the signal Monte Carlo. This method for apply- 
ing the correction is particularly appropriate here 
as the polar-angle, azimuthal and momentum dis- 
tributions and their correlations are very similar 
between the /i pair control sample and the signal 
and background. From these studies, an addi- 
tional correction factor of 0.80 to 0.85 is applied 
to the efficiency as predicted by the Monte Carlo. 
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4. Backgrounds 

There are two broad categories of background 
that require suppression: the background arising 
from non-r+r^ sources, the most significant be- 
ing e+e^ — > /i+/i~(7), and those coming from 
standard r decays. In the latter category, the 
T —^ fih'h'j form an irreducible background. 

The T backgrounds are supressed by requiring 
the signal 7 to be energetic, > iOMeV, and 
by imposing tight requirements on the /i-PID of 
the track. 

The non-T background is reduced by requiring 
there to be a non-/z tag on the non-signal hemi- 
sphere and by removing events with measured 
momentum and calorimeter energy characteris- 
tic of events and hemispheres with little or no 
undetected energy. These events are particularly 
problematic as they naturally possess a value of 
AE that overlaps the signal. 

After all the selection criteria have been ap- 
plied, before applying the signal-box cut, there 
remain background sources in the Grand Side- 
Band region that, from the Monte Carlo simula- 
tion, consist of T ^ fivu (85.9%), t ^ ttv and 
T^Kv (10.6%) and t ^ pv (3.5%). The non-r 
background is expected to be very small in the 
final sample. 

The background is estimated by fitting the 
side-bands in itiec on a sample of data selected to 
have — 2tTA_E < AE < Sctab- As seen in Figure |[ 
the distribution of the background is reasonably 
uniform in the region of the final selection, which 
enables one to use a simple linear interpolation to 
estimate the density of background one expects in 
the signal-box region. A geometrical correction is 
applied to estimate the background contained in 
the elliptical signal-box. 

A low statistics check of the method is made by 
applying it to the r+r" Monte Carlo sample. In 
this case, the luminosity scaled number of events 
predicted from the side-bands is 9.1 ± 2.1 events, 
in good agreement with the 5.5 ± 3.2 events ob- 
served. This estimate of 5. 5 ±3. 2 is the only back- 
ground estimate based purely on counting Monte 
Carlo events. 

The data side-band measurements yield esti- 
mates of 7.8 ±1.4 events. 



5. Systematics 

Using the r+r^ Monte Carlo, the above selec- 
tions yield an absolute efficiency for the signal 
of 6.7 ± 0.1{stat)%. After the final selection of 
requiring the events fall within the ttiec - 
signal-box, the efficiency is 5.2 ± Q.l{MCstat) ± 
Q.b{sys)%. This efficency has a number of sys- 
tematic uncertainties associated with it including 
those arising from: 

1. the trigger efficiency 

2. the tracking reconstruction efficiency 

3. the neutral cluster reconstruction efficiency 

4. the bacfcground filter and skim selection effi- 
ciency 

5. electron opposite hemisphere requirements 
(electron tag) 

6. mopp requirements ('p' tag) 

7. ^-PID requirements 

8. the photon energy scale and resolution 

9. the photon direction reconstruction, scale and 
resolution 

10. the track momentum scale and resolution 

11. the track momentum direction scale and reso- 
lution 

12. the beam energy scale 

13. the beam energy spread 

Evaluation of the efficiency done using the events 
in the Grand Side-Band, as these have character- 
istics which are very similar to those of the signal. 
The good agreement between the data and Monte 
Carlo normalization for the background sources, 
evident in Table 1, and the shapes of, for exam- 
ple the /i momentum, indicates that the efficien- 
cies are well understood. The effects of systematic 
items (l)-(7) are incorporated into the ratio of ob- 
served to expected events: 1.022 ± 0.069(stat) ± 
0.025(norTO). This yields an estimate of 7.3% on 
the systematic error associated with items (l)-(7). 

A global estimate of tracking and calorimetry 
errors is provided by shifting the signal-box posi- 
tion and by modifying the size of the box accord- 
ing to the uncertainties in the resolution. Half 
of the observed changes are used to estimate the 
systematic errors. The systematic errors assessed 
in this manner are presented in Table H. 
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Systematic Influence 
on the Signal Efficiency 


Relative Error 

(%) 


Effects (l)-(7) 


±7.3 


Track and Ecal Resolution: 




/S.E scale 


±0.8 


/S.E resolution 


±3.4 


niEC scale 


±0.3 


mEC resolution 


±0.6 


Ecal scale 


±3.3 


Momentum scale 


negligible 


beam energy spread 


±0.3 


Total 


±8.8 



Table 1 

Systematic Errors associated with the signal ef- 
ficiency. 
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Monte Carlo simulations revealed that the se- 
lection efficiency is insensitive to the uncertainties 
in the beam energy spread and the associated sys- 
tematic error is negligible. 

6. Limit 

From the observed number of events a 90% up- 
per limit is set on B(t ^ fi'y) with the system- 
atic errors included as suggested in using the 
technique of Q. 

A demonstration with the two sidebands can be 
made: the lower sideband had a predicted back- 
ground rate of 6.1 ± 2.2 events, whereas 6 are 
observed. 

When the signal-box was unblinded, we observe 
13 events, shown in Figure I for the mEC projec- 
tion. Of these, 9 had electrons identifed in the 
opposite hemisphere; 5 were h > 7r°-like and one 
was identified as both an electron and a h > 7t°. 
These tagging ratios are consistent with the ex- 
pectations of T decays. The distributions of the 
niEc and AE are shown in Figure ^. 

The background was estimated to be 7.8 ± 1.4. 
This yields a limit of 11.5 event upper limit on 
the number of signal events @90% CL when the 
systematic errors are included. The 11.5 event 
upper limit translates into a limit: B(r — * /17) < 
2.0 X 10-^@90%CL. The probability of a back- 
ground of 7. 8 ±1.4 events up to 13 observed events 
in the absence of a signal is 7.6% if one includes 



Figure 5. Distribution of mEc for the data 
and T+r^ Monte Carlo simulation for the r — > 
/X7 simulated signal for those events within the 
— 2ctab < AE < 3(TAE region. This is after all 
cuts but that on the signal-box. This plot is made 
after unblinding. 



the systematic errors. 

7. Conclusion 

The 1999-2001 BABAR data has been studied 
in a search for the forbidden decay r — > iij. These 
studies reveal that a search with an efficiency of 
5.2 ± 0.1{MCstat) ± 0.5(sys)% and an expected 
background rate of 7.8 ±1.4 results in 13 events 
being observed. This leads to a preliminary limit 
of B(t ^ ^7) < 2.0 ± 10"6@90%CL. 
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